Skip to content

(*)Use por_face_area in zonal_face_thickness#53

Merged
marshallward merged 2 commits into
NOAA-GFDL:dev/gfdlfrom
Hallberg-NOAA:porous_face_thickness
Dec 29, 2021
Merged

(*)Use por_face_area in zonal_face_thickness#53
marshallward merged 2 commits into
NOAA-GFDL:dev/gfdlfrom
Hallberg-NOAA:porous_face_thickness

Conversation

@Hallberg-NOAA
Copy link
Copy Markdown
Member

Use the por_face_area[UV] in the effective thickness calculations in
zonal_face_thickness and merid_face_thickness, so that they are more consistent
with their use elsewhere in the code for the relative weights in calculating the
barotropic accelerations. Because these por_face_area arrays are still 1 in all
test cases, the answers are unchanged in any test cases from before a few weeks
ago, but there could be answer changes in cases that are using the very recently
added capability (in PR #3) to set fractional face areas. This change was
discussed with Sam Ditkovsky, and agreed that there is no reason to keep the
ability to recover the previous answers in any cases that use the recently added
partial face width option.

This commit also expanded the comments describing the h_u and h_v arguments to
btcalc(), zonal_face_thickness(), and merid_face_thickness() routines, the
diag_h[uv] elements of the accel_diag_ptrs type and the h_u and h_v elements of
the BT_cont_type.

All answers and output are bitwise identical in the MOM6-examples test suite
and TC tests, but answer changes are possible in cases using a very recently
added code option.

  Use the por_face_area[UV] in the effective thickness calculations in
zonal_face_thickness and merid_face_thickness, so that they are more consistent
with their use elsewhere in the code for the relative weights in calculating the
barotropic accelerations.  Because these por_face_area arrays are still 1 in all
test cases, the answers are unchanged in any test cases from before a few weeks
ago, but there could be answer changes in cases that are using the very recently
added capability (in PR #3) to set fractional face areas.  This change was
discussed with Sam Ditkovsky, and agreed that there is no reason to keep the
ability to recover the previous answers in any cases that use the recently added
partial face width option.

  This commit also expanded the comments describing the h_u and h_v arguments to
btcalc(), zonal_face_thickness(), and merid_face_thickness() routines, the
diag_h[uv] elements of the accel_diag_ptrs type and the h_u and h_v elements of
the BT_cont_type.

  All answers and output are bitwise identical in the MOM6-examples test suite
and TC tests, but answer changes are possible in cases using a very recently
added code option.
@codecov
Copy link
Copy Markdown

codecov Bot commented Dec 23, 2021

Codecov Report

Merging #53 (2469107) into dev/gfdl (1028ee0) will not change coverage.
The diff coverage is 0.00%.

Impacted file tree graph

@@            Coverage Diff            @@
##           dev/gfdl      #53   +/-   ##
=========================================
  Coverage     28.91%   28.91%           
=========================================
  Files           242      242           
  Lines         71252    71252           
=========================================
  Hits          20605    20605           
  Misses        50647    50647           
Impacted Files Coverage Δ
src/core/MOM_barotropic.F90 37.70% <0.00%> (ø)
src/core/MOM_continuity_PPM.F90 39.46% <0.00%> (ø)
src/core/MOM_variables.F90 28.57% <ø> (ø)

Continue to review full report at Codecov.

Legend - Click here to learn more
Δ = absolute <relative> (impact), ø = not affected, ? = missing data
Powered by Codecov. Last update 1028ee0...2469107. Read the comment docs.

@marshallward
Copy link
Copy Markdown
Member

Gaea regression: https://gitlab.gfdl.noaa.gov/ogrp/MOM6/-/pipelines/14450 ✔️

@marshallward marshallward merged commit 2b2214d into NOAA-GFDL:dev/gfdl Dec 29, 2021
@Hallberg-NOAA Hallberg-NOAA deleted the porous_face_thickness branch January 3, 2022 13:33
@Hallberg-NOAA Hallberg-NOAA added answer-changing A change in results (actual or potential) bug Something isn't working labels Feb 1, 2022
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

answer-changing A change in results (actual or potential) bug Something isn't working

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants